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Abstract 

We consider two sequential models of deposition and aggregation for particles. The 
first model (No Diffusion) simulates surface diffusion through a deterministic cap- 
ture area, while the second (Sequential Diffusion) allows the atoms to diffuse up to i 
steps. Therefore the second model incorporates more fluctuations than the first, but 
still less than usual (Full Diffusion) models of deposition and diffusion on a crys- 
tal surface. We study the time dependence of the average densities of atoms and 
islands and the island size distribution. The Sequential Diffusion model displays a 
nontrivial steady-state regime where the island density increases and the island size 
distribution obeys scaling, much in the same way as the standard Full Diffusion 
model for epitaxial growth. Our results also allow to gain insight into the role of 
different types of fluctuations. 
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1 Introduction 



The diverse morphologies that occur as result of the growth of epitaxial thin 
films are not only of great interest from the point of view of technological 
applications, but also a fascinating subject for basic science. How the inter- 
play between simple microscopic processes (deposition, diffusion, aggregation) 
gives rise to complex patterns on larger scale is a nontrivial question that has 
attracted a lot of interest, in particular since the spectacular refinement of ex- 
perimental techniques has provided an unprecedented control over real growth 
processes [1]. 

The growth of the first monolayer has a special importance: islands form a 
template over which further growth proceeds, and its morphology may have 
a strong infiuence on the nucleation and growth of subsequent layers. For 
this reason submonolaycr growth has been intensely investigated [2] and much 
progress has been done in the comprehension of the different stages in the 
growth process. 

Prom the theoretical point of view, the preferred approach has been the use 
of rate equations [3,4]. They allow to understand well the behavior of mean 
quantities, like the scaling of island density with respect to time, flux intensity 
and temperature [5,6]. However the application of rate equations to charac- 
terize more in detail the morphology of the growing layer, by computing the 
Island Size Distribution (ISD), has turned out to be much more difficult, even 
in the simplest case of irreversible growth. 

Kinetic Monte Carlo (KMC) numerical simulations support the conclusion 
that, instead of depending separately on the coverage 6 and the ratio between 
the diffusion constant of adatoms D and the incoming flux F, the ISD col- 
lapses, in the steady-state precoalescence regime, on a unique function when 
sizes are scaled by their mean value s.^^{9^ D / F). The naive application of rate 
equations with constant capture numbers [5,6] yields unsatisfactory results: 
the ISD scales with s^v, but with a scaling function that diverges for a flnite 
value of its argument and is zero beyond it, at odds with the smooth func- 
tion found in KMC simulations. The origin of the failure has been traced to 
the mean-fleld assumption that the typical environment of an island is inde- 
pendent of its size. Instead, larger islands tend to have larger capture zones 
and as a consequence, larger capture numbers. This has led to very sophisti- 
cated improvements of the original theory in order to overcome this difficulty 
and provide a correct description of the evolution of the system. Bartelt and 
Evans [7] have related the form of the scaling function for the ISD to a generic 
scaling form for the capture numbers. This last quantity is then determined 
via the solution of an evolution equation for the joint probability distribution 
of island sizes and capture zones, under plausible assumptions and simplifi- 
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cations [8,9,10,11]. Other approaches include a general analysis in the limit 
of infinitely fast adatom diffusion [12] and a self-consistent coupling of rate 
equations for island size with equations for capture zones [13]. 

Very recently some doubts have been cast over this conceptual framework and 
in particular about the very existence of scaling for the ISD in the steady- 
state. Authors of Ref. [14] have claimed that scaling of the ISD with respect 
to coverage does not rigorously occur for point island models, while it does 
for extended island models. The ultimate reason for this difference is ascribed 
to the fact that for point islands nucleation continues throughout the steady- 
state, while it is almost completely suppressed for extended islands. According 
to Ref. [14] continuous nucleation would lead asymptotically to a configuration 
where the size of the capture zones is delta-like distributed and this would 
imply the absence of scaling. 

In any case, it is evident that considerable progress has been achieved via the 
refinement of the rate equation approach. In the present work, we take a differ- 
ent path. Instead of trying to refine further the application of the rate equation 
approach to the standard model for submonolayer growth, we investigate, by 
means of suitably modified models, what is the effect of the different types of 
fluctuations that enter submonolayer growth. Work along this line has been 
previously performed in Ref. [15], where it is shown that the ISD obtained 
in KMC simulations is well reproduced by a deterministic continuum model, 
the level-set method [16], with the sole additional ingredient that the spatial 
positioning of new islands occurs stochastically with probability depending on 
the local value of the adatom diffusion held. 

In our paper, we stiidy two simplified atomistic models that switch off some 
sources of fluctuations: this allows, on the one hand, to investigate the precise 
role of each of them, and on the other hand, to understand the origin of the 
scaling property of the ISD. Both models are sequential, in the sense that the 
deposition of an adatom takes place only when the one deposited before does 
not move any more. In the flrst case [No Diffusion model (ND)], a deposited 
adatom remains where it is deposited, unless it is within a flxed distance 
from an island that incorporates it. In the Sequential Diffusion (SD) model, 
an adatom performs a random walk immediately after deposition. If it is not 
incorporated after a flxed number of steps it stops forever, becoming the seed of 
a new island. Both models have already been used in the context of multilayer 
growth [17,18] and the flrst (No Diffusion) model has been recently studied 
for submonolayer growth in one-dimension [19]. 

The two models are in some sense successive approximations of the standard 
(point-island) Full Diffusion (FD) model [5] that is thought to include all 
relevant physical ingredients for irreversible submonolayer growth in the pre- 
coalescence regime. In such a model several mechanisms give rise to spatial 
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and temporal fluctuations: the deposition flux, the diffusion of adatoms, the 
nucleation of new islands. The two models considered in this paper shut off 
some of these noise sources. Temporal fluctuations in the flux are disregarded 
because both models are sequential. The only source of fluctuations in the ND 
model is provided by the spatial fluctuations of the flux; all the rest of the 
dynamics is deterministic. The SD model has an additional type of spatial 
noise, associated to the diffusive motion of adatoms. The main difference with 
respect to the standard model is in the nucleation process, that requires in 
practice only one particle instead of two. 

In the following we study the temporal evolution of mean adatom and island 
densities and the scaling properties of the ISD. It turns out that the behavior 
of the No Diffusion model is far, even from the qualitative point of view, from 
the FD model. Its steady-state is characterized by a complete lack of nucle- 
ation events and, as a consequence, by a "trivial" scaling of the ISD. This is 
due to the fully deterministic treatment of adatom diffusion. The Sequential 
Diffusion model, instead, captures most of the physics of FD model, exhibit- 
ing a steady-state with continuous nucleation and a nontrivial scaling of the 
ISD. The comparison between the two models makes apparent that scaling 
of the ISD can have two distinct origins: either the complete suppression of 
island nucleation in the steady-state or the persistence of nucleation events 
nontrivially coupled to the evolution of capture numbers. The shape of the 
ISD scaling function mirrors the different mechanisms at work. 



2 The Models 



We consider two different models for irreversible submonolayer growth, charac- 
terized by an increasing amount of fluctuations. They are point-island models, 
i.e. the spatial extent of islands is always one lattice unit, island size being 
simply a counter attached to them [5]. If Np particles have been deposited on 
a square substrate of L x L sites, the ratio 9 = Np/L"^ defines the coverage 
9. Point-island models are expected to describe well the behavior of extended 
island models, provided the total coverage is small, ^ <S 1, so that islands do 
not coalesce. 

We call the first model "No Diffusion (ND) model", defined as follows. At 
each time step we randomly choose a lattice site and check whether some 
atoms or islands exist within a distance If they do, the newly deposited 
adatom sticks to the closest atom or island, otherwise it sticks permanently 
to the deposition site. Since adatoms do not diffuse, the distinction between 
adatoms and islands is only formal in this model: adatoms are in all respects 
islands of size one. We keep the distinction for comparison with the standard 
Full Diffusion model. Random deposition is the only source of noise in the ND 
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Fig. 1. Reduced density of adatoms hi (upper lines for small p) and islands his 
(lower lines for small p) for the No Diffusion model with i = 5 (solid), i = 10 
(dashed), £ = lA (dotted), showing that only p matters. Data are for a system of 
size L = 100. 



model. 



The length i is the only length scale in the problem (apart from the lat- 
tice size). Hence all physical quantities are the same provided all lengths are 
rescalcd by £ (as long as £ ^ 1, so that discreteness of the lattice can be 
neglected). Fundamental quantities are the densities, defined as = N^/L"^, 
where is the number of objects of type x (adatoms, islands) present in 
the system. The properly rescaled densities are reduced densities, defined as 
hx = nNx/{L/£y = crix, where c = vrf^ is the capture area surrounding a 
particle. The physical meaning of the reduced density is the number of 
particles of type x per capture area. Also the coverage depends on a length 
(L). The properly rescaled coverage is p = 7rNp/{L/£)'^ = cd, i.e. the average 
number of particles deposited per capture area. The fact that the evolution of 
the system depends on 9 and £ only through the parameter p = cd is verified 
numerically for the reduced densities in Fig. 1. As discussed in Sec. 4.2.1, for 
some properties this model is equivalent to the Random Sequential Adsorption 
model [20]. 

The second model considered is the Sequential Deposition (SD) model. We 
choose at random a lattice site and put there a new adatom. This adatom 
diffuses randomly for at most £ steps. If it happens to go on an occupied site, 
it sticks to the preexisting island (or adatom). If the adatom does not meet 
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Fig. 2. Reduced density of adatoms ni (upper lines for small p) and islands his 
(lower lines for small p) for the Sequential Diffusion model with i = 20 (solid), 
£ = 40 (dashed), £ = 80 (dotted), showing that only p matters. Data are for a 
system of size L = 200. The value of p is computed using Eq. (1). 

any island or adatom during its whole random walk, it stops forever on the 
last site becoming an "island" of size 1. Also in this case we will continue to 
call atoms adatoms, but once again the distinction between such adatoms and 
islands is formal. 



Also in the SD model all quantities do not depend separately on 9 and £, 
rather on the combination p = c{i)6. The definition of the capture area c{i) is 
in this case slightly different from the previous one: it has to be intended as a 
probabilistic concept, the region surrounding an island such that a deposited 
adatom is likely to be incorporated by such an island, before ending its walk. 
Hence the capture area is given by the number A^dis(^) of distinct sites visited 
by a walker during i steps: c = N^i^{i). In two dimensions it is well known [21] 
that, for asymptotically large £ 



log(^) ■ 



The numerical proof that p is the only relevant parameter for the SD model is 
provided in Fig. 2. The sources of noise in this model are the random deposition 
and the diffusion process of adatoms. 

The Sequential Diffusion model bears some resemblance to another model used 
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for describing epitaxial growth in the presence of impurities or surfactants, 
the spontaneous nucleation (SN) model, also denoted by the size i = of 
the critical nucleus [22,23]. The prescription for the SN model is that at each 
diffusive step adatoms have a finite probability to stop, nucleating in this way 
a new island. Intuition suggests that if the probability to stop is such that 
the mean number of steps taken by the adatom is £, the model with i = 
might behave in a fashion very similar to the SD model. We will discuss below 
whether this idea turns out to be correct. 

Since in both the ND and the SD models the only parameter is p = c^, it is 

important to remark that large p docs not imply necessarily large coverage 6. 
With sufficiently large i, the asymptotic stage of these models (see below) can 
occur even for 6* -C 1, i.e. when the assumption of point islands is a reasonable 
approximation for extended islands. 

The SD model incorporates more fluctuations than the ND model, but it still 
lacks two (related) features of the models usually studied in KMC simulations: 
the fact that adatoms are continuously diffusing and the parallel, not sequen- 
tial, character of the deposition process. If we take both into accounts, we 
have the standard model for irreversible growth for point islands [5], that we 
call Full Diffusion (FD) model. A flux F of adatoms is deposited at random 
on the sites of a lattice. Once deposited, each adatom diffuses (with diffu- 
sion constant D) until it meets a preexisting island (thus being incorporated) 
or another adatom (nucleating a new island of size 2). This model has been 
thoroughly studied [5]. We will not simulate it but we will often refer to its 
well known behavior. It is not sequential and it contains all sources of noise: 
temporal and spatial randomness in the deposition process plus the whole 
fluctuations of adatom diffusion and island nucleation. 



3 Behavior of average densities 

The simplest quantities that describe the process of submonolayer growth are 
the reduced density of adatoms fii and of islands fiis. In Figure 3 we report 
the behavior of the average adatom and island densities as a function of the 
rescaled temporal variable p. 

In both cases a crossover at p ~ 1 separates two different regimes. For small 
p nucleation dominates, the density of atoms increases linearly, and the den- 
sity of islands increases quadratically, as predicted by rate equations (see next 
section). At large p the island density cannot decrease, because coalescence of 
islands is not possible in point-island models. For the ND model. His reaches 
asymptotically a saturation value n^. Instead for the SD model the island 
density grows logarithmically before attaining a constant value when the dis- 
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Fig. 3. Reduced density of adatoms hi and islands his for the ND model (a) and the 
SD model (b). Solid lines are the solutions (3) of the rate equations. The decreasing 
dashed lines are power-law decays as p"^/^ (a) and p^^ (b). The increasing dashed 
lines are: (a) a fit to his = V^oo — Ci/-y/p, with Ri 0.45 and ci 2.2 (see 
Sec. 4.2.1); (b) the analytic curve his = (l/a)log(l + ctp), with a ~ 2.93 (see 
Sec. 4.2.2). Data are for £ = 10 (a) and £ = 100 (b). 
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creteness of the lattice comes into play ^ . The density of adatoms decreases 
for large p as for the ND model and as for the SD case. 

These patterns of behavior can be qualitatively understood as follows. At the 
beginning, new atoms are typically deposited in empty regions so that their 
density grows linearly with the coverage. Only occasionally a new particle 
lands in the capture area of a preexisting adatom, giving rise to a new island. 
The average distance among atoms or islands is very large. With increasing 
p, atoms and island densities increase and when p becomes larger than 1 
their capture areas start to overlap. Prom then on the growth processes of 
distinct islands are no more independent and the difference in the models 
results in a different evolution. At large p, almost all neighbouring islands 
have overlapping capture areas but they do not cover the whole lattice. There 
are still some regions such that all their sites are outside the capture area of 
any particle, so that the new adatoms deposited there do not get incorporated. 
Densities of atoms and islands change because these active regions are filled 
by adatoms and adatoms are then transformed into islands. While for the 
ND model capture areas are deterministic, so that once they cover the whole 
lattice no more nucleations are possible, for the SD model it is always possible 
that an adatom deposited close to an island does not in fact meet the island 
during its diffusive motion. Such adatom will stop at the end of its walk and 
give rise to a new island. This difference explains why the island density tends 
to a constant for the ND model, while it keeps growing (though slowly) in the 
SD case. In the latter case the growth of fiis ends only when the island density 
is so high (of the order of 1) that new particles can only be incorporated in 
existing islands. 

These results must be compared to the behavior of the analogous quantities 
for the FD model. When adatom diffusion is fully operative, the density of is- 
lands grows as 9^ during the initial regime, while the adatom density increases 
linearly. In the ensuing steady-state regime instead, grows as (FO/DY^^. 
Simultaneously the density of adatoms decreases as (FOD^)^"^^^ [5]. From this 
comparison it turns out that, from the quantitative point of view, both mod- 
els with reduced fluctuations fail to reproduce, at least for what concerns the 
exponents, the phenomenology of the FD model. For example, the temporal 
exponent for the growth of riis during the nucleation stage is 2, instead of 
3. This is a direct consequence of the fact that in the FD model two mobile 
adatoms must meet to nucleate an island, while in the present models nu- 
cleation involves a single adatom, that, if not immediately incorporated after 



The transition to the final asymptotic stage for the SD model with constant 
island density occurs only for values of p ^ c. In such a limit, corresponding to the 
condition 9^1 for which point island models become unphysical, the adatom and 
island densities approach exponentially and a value of the order of 1, respectively. 
In our simulations we have never reached such extremely high values of p. 
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deposition, becomes immobile, in this way seeding a new island. However, from 
a more qualitative point of view, the SD model performs rather well, since it 
exhibits regimes which have the same nature of the more realistic FD model. 
In particular, a steady-state regime during which islands are nucleated is non- 
trivial. This means that, after a sufficiently long time interval, most islands 
have been nucleated during this regime and the nucleations occurred during 
the initial stage are negligible. This is what happens also in the FD model and 
will be shown to have profound consequences for the scaling of the Island Size 
Distribution. 

In the stationary-state the behavior of the density of islands is then the direct 
consequence of the amount of spatial fluctuations present in the nucleation 
process. The failure of the ND model in reproducing such a regime indicates 
that the existence of an area around particles where new nucleation events are 
deterministically forbidden is a too drastic simplification. 

In the two models considered here, the quantity I determines the value of the 
densities in the steady-state. Hence it plays a role analogous to D/F in the 
FD model. Another difference deserves to be remarked. In the Full Diffusion 
model the morphology of the system depends on two parameters, 9 and D/F, 
that can be changed independently. In the ND and SD models instead, 9 and 
I are combined in the only relevant parameter p. This means that the scaling 
property of the ISD (see below) holds with respect to both parameters or to 
none. 



4 Rate equations 



4.-1 Nucleation stage p -C 1 



As discussed above, for both the ND and the SD model particles are effectively 
"non-interacting" in the limit of small p, because their mean distance is so large 
that typically capture areas do not overlap. In this limit it is easy to write 
down the following mean-field rate equations: 



^ = 1 - (n,, + 2n0 (2) 
dp 

dfijs _ . 

dp 

with initial conditions ni(0) = ^^^(O) = 0. Notice that the rate equations are 
the same for the two models, the only difference being the definition of p. The 
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fact that the rate of island nucleation is proportional to hi reflects the fact 
that only one adatom is necessary for it. 

The solutions are 



^is = 1 - (p + 1) exp(-p). 

These expressions, that make sense for any p, are reported in Fig. 3 as full 
lines and compared with numerical results (symbols). The comparison is very 
good up to p ~ 1, i.e. during the initial nucleation stage. For larger values 
the overlap of capture areas becomes relevant and the above approximation is 
no longer valid, as expected. The equations for the evolution of the density of 
islands of any size s > 2 are 



They can be solved recursively, showing that at small coverage the size distri- 
bution is Poissonian 



It is worth noting that the above rate equations, and therefore their solu- 
tions, do not depend on the space dimension. Changing the dimension simply 
changes the definition of the capture area c. 



4-2 Steady-state p ^ 1 

For large values of p the form of rate equations must be modified to take into 
account the overlap of capture areas. Rate equations for a generic model of 
irreversible growth are 




(3) 




(5) 



dhi 



(6) 



dp 
dhs 
dp 



^agg(s - 1) 



^agg(s) 



S>1 



S > 1 
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where R^gg{s) is the reduced aggregation rate for islands of size s. By summing 
over s, the equation for the total density of islands fiis is obtained 

^ = ^..(1). (7) 



The specification of the aggregation rates R^gg{s) fully determines the set of 
equations. In analogy with the procedure for the Full Diffusion model [7] we 
define the capture numbers cr^ by the relation 

-Ragg(s) = CTsfls- (8) 



When all capture numbers are equal to one, wc recover Eq. (2) valid during 
the nonintcracting nucleation stage. In Section 5, devoted to the study of 
the Island Size Distribution, we will discuss the appropriate expressions for 
capture numbers as in the regime for large p and make use of Eq. (6). Here we 
study the asymptotic behavior of the average densities by means of approaches 
tailored on the specific model. 



4-2.1 No Diffusion model 

When p is large, almost all lattice sites are at a distance d < £ from the 
closest island, i.e. they are part of a capture area. However, there are some 
domains whose sites are aX d > i from every island or adatom: we denote as Iq 
the number of such active regions. Similarly, we call /i the number of regions 
corresponding to the capture zones of adatoms. Finally, let us call the 
average size of capture zones for p ^ oo, and — A^/c the corresponding 
reduced quantity. 

Let now Np be the number of particles deposited and Ay {9) the mean size of 
the active regions. At very high p we can assume that the deposition of an 
adatom in an active region (occurring with probability lQAy{6)/ L'^) results in 
the complete cancellation of such active region, so that Jq is reduced by one. 
If instead an adatom is deposited in a region close to a preexisting adatom 
(and this happens with probability Aooh/L'^) a new island will be nucleated 
and /i — > /i — 1. Hence the following equations hold 



dip _ Aylp 
dNp ~ L2 
dh ^ Aylp A^ 
dNp~ L2 L2 
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s 

Fig. 4. Distribution of the size of active regions at p ~ 245. Symbols: numerical 
results for the ND model with ^ = 10. Line: analytical results according to Eq. (11). 

The crucial point is the determination of the form of Ay for p ^ 1. For large 
p active areas can be considered to be independent. The evolution equation 
for the number N{s,p) of active areas of size s is dN{s,p)/dp = —sN{s,p)/c 
so that 

N{s,p) ^ N{s,po)exp[-s{p-po)/c] (10) 



where Pq is the value of p (of the order of 1) after which active areas can be 
considered as independent. The form of N{s,Pq) depends on the dynamics 
during the nucleation stage for small p. For its determination it is useful to 
take advantage of the similarity between the No Diffusion model and the so 
called Random Sequential Adsorption (RSA) model, that is well studied in 
the literature [20]. Our model can be mapped onto RSA for disks of radius 
r = £/2. There are some differences in the temporal properties of the two 
models, but for what concerns spatial properties the mapping is exact. We 
can use the result for the size distribution [24] 

Ar(s,po)oc(s-1)-V2 (11) 



(see Fig. 4 for a numerical check), so that 

jrsNis,p)ds 1 
Ay^is)^ /^oo^(^^^)^^ -1 + ^- (12) 
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The first term appearing in the r.h.s. of Eq. (12) reflects the discreteness of 
the lattice, the obvious fact that no region can be smaller than 1. When active 
regions are on average much larger than 1 site, the first term can be neglected. 
This occurs for ^ ^ 1 (even if p » 1, because I is large). Inserting At/ 1/(2^) 
into Eqs. (9) yields 



^ = (13) 
dp 2p ^ ^ 

dp 2p °° 

As expected everything depends only on p = c{t)9. The solutions are 



(14) 



Hp) 



2^ fe=i 



{2k - 1) 



-1 pk+\ 



If instead ^ ^ 1, ^y{0) Ri 1: all remaining active regions have very small size 
and lattice discreteness dominates. One finds that ni decays exponentially to 
zero and in the same way Uis goes to its asymptotic finite limit. In this limit 
observables depend separately on £ and 0. 

For p — > oo Equations (14) imply 



hi{p) (X h 



Cl 



2A^I^'^ 



(15) 



The power- law decay with exponent 3/2 for the adatom density is in agreement 
with simulations (see Fig. 3). Since for p — > oo both Iq and Ii vanish, only 
islands remain in the system and their capture zones cover the whole lattice. 
The asymptotic value of their density is then n^s = 1/Aoo with corrections, 
for finite p, proportional to the largest of the two vanishing contributions. 



4-2.2 Sequential Diffusion model 

In the Sequential Diffusion model, a freshly deposited adatom wanders until 
it sticks to a preexisting adatom, to an island or it stops after I random-walk 
steps. Let us denote as pi, pis and po the normalized probabilities for these 
events to occur, respectively. The evolution of the average densities for the 
Sequential Diffusion model can be written as 



14 



dhi 
dp 
dhis 
dp 



(16) 



Notice that since preexisting adatoms and islands are perfectly equivalent for 
the walker, the ratio pi/ Pis must be equal to the ratio of the densities ni/fiis, 
so that 



Pi = (1 - Po) 



Pjs = (l -Po)- 



hi + hi 



Til + His 



(17) 



This implies that the equation for the total density of objects ht = hi+ his is 
simply dht/dp = Pq. This last quantity is nothing else than the survival prob- 
ability after i steps for a random walk in a lattice with a density rit of traps. 
In the initial regime p the spatial distribution of traps is uncorrelated. It 
is well known that [21] for small trap density rit, 

Po^{l-ntf'^-. (18) 



We remind that the number N^^^ of distinct sites visited by a walker during 
i steps is by definition the capture area c. Using Eq. (18) and the condition 
nj -C 1, Eqs. (17) become the rate equations (2) valid for small p, as expected. 
The total density obeys then 

dhi , , , , 

= exp{-ht) (19) 



with solution 

ht = \og{l+p). (20) 



Equation (20) is exact for p 0, but it makes formally sense also for large p 
where it predicts logarithmic growth, in qualitative agreement with numerical 
results. However, a close comparison shows a quantitative mismatch (Fig. 5). 
The origin of this difference is clear: we have assumed an uncorrelated dis- 
tribution of objects. This is appropriate at the beginning of the dynamics, 
but it becomes incorrect when the overlap between capture areas starts to be 
significant. To overcome this difficulty, let us study G{x), the distribution of 
distances between nearest neighbor objects (both islands and adatoms). 
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Fig. 5. Behavior of ht for the SD modeh the square symbols are the results of nu- 
merical simulations for i = 100; the solid line is the theoretical prediction according 
to Eq. (20), the dashed line given by Eq. (24) with a 2.9. 

Let us consider first the distance distribution G„„^(a;) for objects put in the 
lattice without correlations (with probability rit) . In order for an object to be 
at distance x from its nearest neighbour, the i/(x) sites at distance y < x must 
be empty: this happens with probability l — rit for each of them. Moreover, at 
least one of the S{x) sites at distance x must be occupied. Then 

G_(x) oc (1 - ntY^^^ • [1 - (1 - ntY^^^] (21) 



The quantity u^x) can be reasonably approximated with the continuum ex- 
pression i>{x) — TTX^, resulting in a Gaussian contribution to G^^^{x). The 
other factor, equal to 27ix in a continuous approximation, strongly fluctuates 
for the different integers, due to the discreteness of the lattice. In Fig. 6 we 
compare G{x) with G„nc{x) for small (a) and large (b) values of p. While for 
p -C 1 the agreccmcnt is excellent , for large p there is a systematic deviation 
for large distances. 

To understand which part of the distribution G{x) of distances between nearest 
neighbors is relevant for the determination of the probability po that enters 
Equations (17) for the SD model in the late stage, let us consider first for 

^ The agreement is very good but not perfect, even for p ^ 0, because for dis- 
tances smaller than ^/c correlations are always at work and one expects a reduced 
probability G{x) with respect to the fully uncorrelated case. 
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simplicity the ND model in one dimension. In such a case, po is the probability 
that a newly deposited particle is at least at a distance i from all islands in the 
system. Contributions to po are given only by the intervals between particles 
larger than 2i so that 



For the SD model, due to the stochastic nature of diffusive motion one cannot 
write down such a simple exact expression. However, Eq. (22) remains approx- 
imately true provided the capture area 2£ is replaced with the capture area for 
the SD model N^i^{i). In two dimensions additional complications arise, be- 
cause of the angular constraints. Nevertheless it is intuitively clear that also in 
two dimensions the chance for an adatom not to be incorporated depends on 
the existence of regions that are sufficiently far from all particles, and hence 
from the shape of the function G{x) for large x. 

In order to fit the large- ,t decay of G{x) with a form like Eq. (21) one must 
multiply the density by a factor a{p) that is close to 3 for large p (see Fig. 6). 
Qualitatively it is possible to understand why a > 1. The spatial distribution 
of islands obtained with the SD model differs from the spatial distribution of 
the random deposition model in one important respect: the probability to find 
very close islands is strongly suppressed by the diffusion/aggregation process, 
while it is significant in the random deposition model. These neighboring is- 
lands do not increase very much the capture probability, because they behave 
as a single trap. In simple words, for the same value of p, the islands of the 
SD model are more effective to trap diffusing atoms than the islands of the 
random deposition model: therefore, the latter model requires a higher density 
of islands (i.e. a higher p) in order to get the same trapping effect. 

If now we take a to be constant, Eq. (19) becomes 

— = exp[-ant) (23) 
dp 

whose solution is 



Such expression shows a much better agreement with numerical results, as one 
can see in Fig. 5. When p becomes large enough, the density of adatoms is 
much smaller than that of islands. Hence 




(22) 



2e 



fit^ - ■ log(l + ap). 



(24) 



a 



Uis ^ht^ -log(l + ap) 



a 



(25) 
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• log[G(x)] 
□ log[G (x)] 

unc^ 

« log[G (x)](a = 2.93) 



(b) 



_L 



_L 







3 4 
X 



Fig. 6. (a) Comparison of G{x) measured in simulations for the SD model (filled 
circles) with the form predicted assuming lack of correlations G^n^{x) [Eq. (21)] 
(empty squares) for p = 0.02. (b) Comparison of G{x) measured in simulations for 
the SD model (filled circles) with the form predicted assuming lack of correlations 
Gunc{x) [Eq. (21)] (empty squares) and G^^^x) with a ^ 2.9 for p = 60 (empty 
diamonds). 
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and, using Eq. (17) and the fact that p^^ ~ 1, 

duis rii 

^ pic^ — (26) 



dp n, 



IS 



from which we obtain 

_ dhis 1 log(l + ap) 

m ~ riis-T— ^ ^— -■ (27) 

dp a 1 + ap 



In this way, we have recovered that, for p 3> 1, ni ~ p ^ (apart from a 
logarithmic correction) and Uis is logarithmic in p (see Fig. 3). 



5 Island Size Distribution 



During the initial nucleation stage the island size distribution Us is given by 
Eq. (5), reflecting the lack of correlations in the growth of different islands. 
In the foflowing steady-state regime the competition between different islands 
strongly perturbs this form. In the Full Diffusion model it is believed [5] that 
the ISD obeys asymptotically the scaling form 

^^(^) = ;^/[^/^-(^)] (28) 



where 



is the average size of an island. In Figs. 7 and 8 we show that the same scaling 
property holds for both models considered here, although the approach to the 
asymptotic shape is much slower for the Sequential Diffusion model. 

While the scaling property is the same, the scaling function is rather different 
in the two cases and this reflects the totally diverse origin of scaling in the 
two cases. 

For the No Diffusion model, scahng is a trivial consequence of the lack of 

nucleation events in the steady state. When the steady state is reached, the 
whole lattice is partitioned into capture zones. In the present case, owing to 
the deterministic nature of the incorporation process, capture zones exactly 
coincide with Voronoi Cells. During the steady-state, all deposited atoms are 
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Fig. 7. Collapse plot of the Island size Distribution for the ND model. Data are for 
£ = 10. 
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Fig. 8. Collapse plot of the Island size Distribution for the SD model. Data are for 
I = 100. 
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immediately incorporated by the island owning the capture zone wherein the 
atom is deposited. Hence the growth rate of an island is (apart from fluctu- 
ations in the deposition flux) perfectly proportional to the size of its capture 
zone. If we wait long enough the island size distribution will then become equal 
to the size distribution of capture zones. In this sense scaling is trivial in the 
ND model. What is highly nontrivial is the shape of the scaling function which 
is the outcome of the complicated evolution during the nucleation stage [25]. 

As shown in Fig. 3, in the SD model island nucleation continues during the 
steady-state, and if we wait long enough virtually all islands have been nucle- 
ated during such p ^ 1 regime. It is clear that in this case the origin of the 
ISD scaling property is different and has to do with the dynamics during the 
steady-state. The situation is analogous to what happens in the Full Diffusion 
model. 

More insight into the problem is provided by the application of the approach 
that Bartelt and Evans developed for Full Diffusion models [7] . The basic idea 
is that scaling of the Island Size Distribution is a consequence of a scahng 
relationship valid for the capture numbers cr^ entering the rate equations (6) 




Introducing the scaling ansatz for Ug and ag into the second of Eqs. (6), and 
using the steady-state condition, that in this case is a^^riis = 1, a relation 
between the scahng functions / and C is obtained [7] 



fix) = /(O) exp 



X 

I 



dy 



2uj 



- C'{y) 
C{y) -ujy 



(31) 



where 



d In s^-^ 

dine 



(32) 



The analytical computation of the scaling form C{x) for capture numbers is 
a hard problem that has been approached in many ways [8,9,10,11]. Here we 
limit ourselves to its numerical evaluation, performed in a way analogous to 
Ref. [7]. Results are displayed in Fig. 9. 

It is clear that capture numbers for the SD differ much from the mean-field 
assumption that they are all the same, i.e. independent of s. More importantly, 
the scaling form (30) tends to be fulfilled, but only in the limit of very large 
p. This goes along with the observation that in the steady-state regime Sav 
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Fig. 9. Main: Scaling plot of the capture numbers according to Eq. (30) for the SD 
model. Inset: Scaling plot of the capture numbers for the ND model. 



grows linearly with 6 but with a (sizeable) logarithmic correction, so that uj in 
Eq. (31) tends to 1 with respect to 9 only in the limit 9 oo. In particular, 
a linear fit for s/Sav between 0.8 and 1.6 in Fig. 9 yields the effective values 
u = 0.40 for p = 55, u = 0.70 for p = 550, lu = 0.81 for p = 3500 and 
LU = 0.83 for p = 5500. We conclude that deviations from scaling for the ISD 
displayed in Eq. (8) are a consequence of the logarithmic corrections. However 
the logarithmic growth of rij^, i.e. the fact that nucleations continue during 
the steady-state, is also the very origin of the scaling observed for SD. In the 
absence of nucleations one would have uj = 1 and C{x) = x, as it happens for 
the ND model (see the inset of Fig. 9). If we insert such values in Eq. (31), the 
integrand is indeterminate, confirming that the approach of Ref. [7] cannot be 
applied for ND, since the shape of the scaling function f{x) is not determined 
by the dynamics during the steady-state. 



Finally, let us remark that the different nature of ISD scaling is also reflected 
in the form of the scaling function f{x). In the ND case, since nucleation is 
inhibited there are no small islands, and /(s/s^v) vanishes for small argument. 
Instead in the SD case, small islands are continuously created and this implies 
that /(O) is finite. 
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6 Conclusions 



In this paper we have analyzed, by means of numerical simulations and theo- 
retical arguments, two models for irreversible submonolayer growth that incor- 
porate less fluctuations than the usual Pull Diffusion model. Let us conclude 
with some remarks on the findings. 

As already mentioned in the Introduction, the definition of the SD model seems 
similar to the model with critical size i = 0, also said spontaneous nucleation 
model. Surprisingly, all results indicate that the two models exhibit quite 
different behaviors. In both models the trivial nucleation regime is followed 
by a steady-state during which island nucleation continues, but in the SD, 
island density grows logarithmically, while it grows as a power law for SN [26] . 
The difference is even larger for the Island Size Distribution. At odds with the 
result [23] that the scaling function for spontaneous nucleation is a decreasing 
function of the scaled size, i.e. its maximum is for s/Sav 0, in the SD model 
the scaling function has a shape with a peak for finite s/s^^, much more similar 
to the Full Diffusion model. This rather strong change can stem only from one 
of the two differences between the models: the SD case is sequential (while 
the SN is not); in the SD a particle diffuses for a fixed number of steps before 
stopping, while this number fluctuates (with a well defined average) for SN. 
The first difference is for sure irrelevant: also the SN model becomes sequential 
in the limit D/F ^ oo and its behavior does not change. The key point is 
the second difference and in particular the fact that in the SN there is high 
chance for an adatom to stop and nucleate very early after its deposition. This 
facilitates the nucleation of new islands also in the steady-state and enhances 
the probability that newly formed islands have very small capture zones and 
therefore relatively small sizes. The fact that the walker must diffuse for £ 
steps turns out to be crucial for having a peak for flnite s in the ISD of the 
SD model. 

Both point-island models analyzed here display scaling of the Island Size Dis- 
tribution in the steady-state. This result is in apparent contradiction with the 
recent statement [14] that point-island models cannot obey scaling, because 
of continuous nucleation leading to a delta-like distribution of capture zone 
sizes. According to this argument, scaling could occur only if nucleation of 
new islands is suppressed during the steady-state, so that each island collects 
a number of adatoms proportional to its surrounding capture zone and the 
scaling function is equal to the distribution of capture zone sizes. 

The ND model obeys scaling, despite being a point-island model, but this is not 
in real disagreement with Ref. [14]. The scaling observed is of the same type: 
it simply reflects the proportionality between capture zones and aggregation 
rates. It is associated with a linear growth of the average island 
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function of coverage and to a scaling function f{x) that vanishes for a; ^ 0. 
This may happen for point-islands because the particular deterministic nature 
of capture areas in the ND model fully inhibits nucleations in the steady state. 

As shown in Fig. 8, also the SD model obeys scaling. This is a different type 
of scaling, that occurs just because nucleation events occur during the steady- 
state. This scaling is characterized by a sublinear growth of the island density 
riis and by a finite value for /(O) and it is the same type of scaling that occurs 
for the FD model. 

The origin of the two types of scaling is completely different: while in the ND 
case the scaling function is determined by what occurs during the nucleation 
stage, in the SD any feature of the nucleation stage is completely erased by 
the steady-state dynamics, the only thing that matters. The second type of 
scaling is compatible with point-island models, because the steady-state can 
be made infinitely longer than the nucleation stage so that all memory of the 
previous regime is lost, still keeping the coverage smaller than 1. In the SD 
model this occurs considering l/c{i) -C ^ <S 1. In the FD model this occurs 
when [D/F)'^/'^ ^ ^ <^ 1. It would be interesting to gain a greater insight 
into the SD model, trying to determine analytically the form of the scaling 
function C and formally prove that scaling occurs. 
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